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Experiments of pure tensile fracture in brittle gels reveal a new dynamic oscillatory instability 
whose onset occurs at a critical velocity, Vc = 0.87Cs, where Cs is the shear wave speed. Until Vc 
crack dynamics are well described by linear elastic fracture mechanics (LEFM). These extreme speeds 
are obtained by suppression of the micro-branching instability, which occurs when sample thicknesses 
are made comparable to the minimum micro-branch width. The wavelength of these sinusoidal 
oscillations is independent of the sample dimensions, thereby suggesting that a new intrinsic scale 
exists that is unrelated to LEFM. 

PACS numbers: 46.50.+a, 62.20.Mk, 47.54.+r 



The dynamics of rapidly moving cracks are of both fun- 
damental and practical importance. Driven by the elastic 
energy stored in a stressed body, brittle cracks acceler- 
ate towards their theoretical limiting speed [l[. How- 
ever in tensile fracture, the crack tip becomes unstable 
at less than half that velocity (Rayleigh wave speed - 
Vr), sprouting small side cracks termed micro-branches 
0| • At this point, the single-crack state is replaced by 
a multi-crack one and much of the energy flowing to the 
crack tip, is transformed into additional fracture surface 
(roughness). The density of micro-branches is observed 
to increase with the main crack velocity and as a con- 
sequence, cracks are rarely observed at velocities greater 
than 0.6 - 0.7Vr. 

We describe experiments performed on polyacrylamide 
gels. Fracture dynamics in these materials have been 
shown || to be identical to those of standard brittle 
amorphous materials, with the advantage that the ex- 
tremely slow wave speeds in gels simplify the study of 
crack dynamics. We demonstrate that, when micro- 
branches are suppressed, crack dynamics are well- 
described by the single-crack equation of motion up to 
previously unattainable crack speeds. Surprisingly, we 
find that a new oscillatory instability occurs at a criti- 
cal velocity close to Vr. The existence of this instability 
highlights important unresolved questions regarding the 
path selection of dynamic cracks Furthermore, the 

suppression of micro-branches when the gel thickness ap- 
proaches the process zone scale, raises interesting ques- 
tions regarding both the origin of the micro-branching 
instability and the relevant length scale at which 3D frac- 
ture becomes effectively 2D. 

Our experiments were performed using sheets of brittle 
polyacrylamide gels of typical size (X x Y) 125 x 115mm 2 
and thickness 0.15 < d < 0.35mm where X, Y and Z are, 
respectively, the propagation, loading, and sheet thick- 
ness directions. All experiments reported here were per- 
formed in polyacrylamide gels composed of 13.8% to- 
tal monomer concentration and 2.6% bis-acrylamide as 
cross-linker 0. These gels have a measured shear mod- 
ulus of G = 35.2 ± lAkPa yielding shear and longitu- 



dinal wave speeds of Cs = 5.90 ± 0.15ra/s and Cl = 
11.8±0.3ra/s, respectively. Measurements of Cs and Cl 
using high-speed photography agree with these values to 
within 0.1m/ s. Experiments using other gel composi- 
tions (e.g., 23% total monomer concentration with 10% 
bis-acrylamide whose Cs ~ 17m/ s) exhibit similar qual- 
itative and quantitative behavior. The gels were cast 
for each experiment between flat glass plates. To ensure 
uniformity, the evaporation of the gel's water content was 
limited to below 5%. We compensated for any resulting 
variations in wave speeds by direct measurement of G for 
each experiment. 

The experimental setup and measurement techniques 
are described elsewhere |3 and will be briefly reviewed 
here. The gel sheets were loaded in Mode I (uniaxial 
extension) via constant displacement in the vertical (Y) 
direction. The applied force was monitored through a 
load cell. When a desired stress was reached, a short 
cut was made at the sample's edge, midway between the 
vertical boundaries. This seed crack either immediately 
or after a short creep accelerated to dynamic velocities. 
Throughout each experiment the crack tip's profile and 
immediate location were monitored by a high speed cam- 
era set to a X x Y spatial resolution of 1280 x 64/96 pixels 
(frame rate of 7. 2 /5kHz). At the completion of the frac- 
ture event, the fracture surface (XZ plane) was analyzed 
via an optical microscope and the fracture profile (XY 
plane) was scanned with 2400dpi resolution. 

The dynamics of a typical experiment in thin gels are 
presented in Fig. la. The initial crack propagates along a 
straight trajectory until losing its directional stability to 
crack tip oscillations in the XY plane. Straight-line prop- 
agation (Fig. lb) typical of micro-branch free fracture 
gives way to sinusoidal oscillations at a millimeter scale. 
The fracture surface formed by oscillating cracks is typi- 
cally smooth (mirror-like), with rare instances of short di- 
rected lines of successive micro-branches ( "branch-lines" ) 
|jj . This behavior is in contrast to propagation in thicker 
samples where the crack retains its straight line tra- 
jectory throughout the fracture process (Fig. lc). In 
these latter events, the fracture surface is micro-branch 
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FIG. 1: (a) A sequence of photographs, shown at 0.69 ms 
intervals, of a propagating crack undergoing a transition from 
linear (top 2 pictures) to oscillatory motion. Photographs of 
XY profile (top) and (XZ) fracture surface (bottom) of (b) 
a 0.2 mm thick gel sample where oscillations developed and 
(c) a 2.0 mm thick gel where the crack retained its straight 
line trajectory. In (c) the fracture surface is micro-branch 
dominated whereas in (b) the oscillating crack has a mirror 
surface. Propagation in (a), (b) and (c) was from left to right. 
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FIG. 2: (color online) (a) Oscillations appear only beyond a 
critical velocity of Vc=0.S7Cs- Vc is independent of the crack 
length at the onset, l c , sample size, L, and thickness, d. The 
different symbols represent measurements for 0.15 < d < 0.30 
mm. (inset) The velocity components, V x (continuous line) 
and V y (dotted line), as a function of X, where the dashed line 
marks the onset of oscillations, (b) The measured (squares) 
and predicted (line) relation between l c and E/a 2 in the 
framework of LEFM The line's slope gives the fracture 
energy T(v c ) = 34 ± 2J/m 2 . 



dominated; beyond a critical velocity branch-lines de- 
velop whose density increases with the crack velocity Q . 
Branch-lines disappear upon arrival at the sample sur- 
faces at Z=0 or Z=d. Thus, in thin sheets branch-lines 
are quickly suppressed and sporadic. 

The onset of the oscillations occurs at a critical ve- 
locity of V c = 0.87 ± 0.02Cs (Fig. 2a). The ob- 
served value of the critical velocity is independent of d for 
0.15 < d < 0.3mm. In thicker gels, cracks rarely reached 
Vc due to micro-branching (more on this later). Surpris- 
ingly, no evidence of the Yoffe instability |g, in which a 
crack is expected to deviate from its original straight line 
trajectory at V > 0.6Vr, is observed. 

The rapid acceleration to the extreme crack velocities 
attained (V > 0. 87 'Cs) suggests that waves reflected from 
the system's remote boundaries can not affect crack dy- 
namics, as reflected waves only catch up with the crack 



tip well beyond the critical crack length I = l c . Thus, 
for V < Vc, predictions by linear elastic fracture me- 
chanics (LEFM) of the equation of motion for a sin- 
gle edge crack in a semi-infinite body [l| should hold: 
V = V R {1 - T{v)E/K 2 ) where K=1.12o"VVZ and T is the 
material's fracture energy At V = Vc this predicts 
L to be: 



r(v c )E 



1.12W(1-V C /V R ) 



(1) 



Eq. [l] predicts that l c oc a~ 2 where a is the applied 
stress. This prediction (Fig. 2b) is, indeed, in excellent 
agreement with measurements of l c . The slope of the 
linear fit yields a value of the fracture energy, T(Vc) = 
34 ± 2J/m 2 , which agrees with previous measurements 
fin| . The slight scatter in Fig. 2b is consistent with 
our 5% allowed variation of G. This is the first direct 
experimental verification of the equation of motion of a 
single crack at velocities approaching Vr [TH - 

At Vc the behavior of the crack's two velocity compo- 
nents, V x and V yi changes abruptly. Prior to Vc, V y = 
whereas beyond Vc, V y oscillates at a wavelength, A(x), 
with increasing amplitude. The oscillations in V y are ac- 
companied by a much smaller oscillatory component of 
V x having a wavelength of A/ 2 and a tt/2 relative phase 

shift. (Fig. 2a inset). This suggests that V = ^jv 2 + V 2 

does not, itself, oscillate but increases monotonically, al- 
though we lack the temporal resolution to directly verify 
this. 

V x increases steadily until reaching Vc- Beyond Vc, 
the full velocity, V, generally continues to accelerate (see 
e.g. Fig. 3) to steady state velocities. "V, however, has 
not been observed to surpass Cs UM • Unexpectedly, the 
mean value, over an oscillation period, of V x remains 
nearly constant (even dropping by a few percent) upon 
reaching Vc- This is illustrated in Fig. 3 (inset) where 
measured values of V are identical to values obtained by 
integrating the measured path length of the crack and 
assuming that V x is constant, equal to its mean value for 

v> v c m. 

Let us consider the nature of the oscillations. The oscil- 
lations are generally sinusoidal except for large amplitude 
states at the highest values of applied stress, a. Both 
the oscillation amplitude, A, and wavelength, A, develop 
with crack advance until reaching steady-state values af- 
ter a number of cycles (Fig. 4 insets). The steady-state 
values of both A and A vary with a, as shown in Fig. 
4, and do not appear to be well-defined functions of V. 
At V = Vc, A initiates from a finite value, although 
no hysteresis in the transition to the oscillatory state 
has been observed. A is a sharply decreasing function 
of a. Steady-state wavelengths range from 4-9mm with 
V x /X ~ 0.5 — 1kHz. This frequency scale corresponds to 
the same (~ 1ms) timescale that was previously observed 
for micro-branching activation times in similar gels 0. 
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FIG. 3: (color online) The full crack velocity, V (continuous 

line), is compared to V x (dotted line). Beyond Vc, 1 

until Cs (dashed line) 12], while V x stays nearly 

(inset) A comparison of direct measurements of V 

(full circles) with V, calculated using the spatial tr 

the oscillations (squares) The dashed line denotes Us. 

All velocities are averaged over an oscillation period. 



The origin of these length/frequency scales is currently 
unknown. The time scale is over an order of magnitude 
smaller than any possible resonance with the system's 
boundaries. In fact, these scales are wholly independent 
of any of the physical dimensions of the sample, since gels 
of substantially different widths and dimensions all fall on 
the same steady-state amplitude/wavelength curve (Fig. 
4). Furthermore, as indicated by Fig. 2b, the crack is 
effectively travelling in a semi-infinite medium at the on- 
set of oscillations and, under semi-infinite conditions, the 
crack will have no knowledge of the actual size of the sam- 
ple. We, therefore, surmise that the length scale must 
have its origins in an intrinsic scale, such as the process 
zone, that is not quantitatively accounted for in current 
theoretical frameworks. 

Although the steady-state values of A and A are well- 
defined, their initial transients (Fig. 4 - insets) are ex- 
tremely diverse, with no discernable relation to either 
V or a. It is interesting, however, that these transients 
follow a number of distinct A — A trajectories, suggest- 
ing that there may be a number of meta-stable solution 
branches that eventually evolve to the steady-states pre- 
sented in Fig. 4. 

Why is the oscillatory instability not readily observed 
in thick samples? Basically, in thick samples cracks rarely 
attain Vc, due to the onset of micro-branching. Let us 
momentarily digress and consider the relation of micro- 
branch formation to sample thickness. Recent work in 
gels indicates that the micro-branching instability is a 
first-order phase transition in which micro-branch for- 
mation is triggered by noise, above a critical velocity Q. 
Micro-branches in gels have a minimal width Q, AZ, 
whose scale (here ~ 40/im) is approximately that of the 
process zone. The system's quenched noise level may 
then be related to d/AZ. This may explain why the 



FIG. 4: (color online) Steady-state amplitudes, A, (a) and 
wavelengths, A, (b) as a function of the applied stress, a. 
(insets) Transient time evolution of A and A in two typical 
experiments. t=0 defines the onset of oscillations. Sample 
dimensions: (X x Y) ~ 125 x 115mm 2 (■), 130 x 155mm 2 
(A) and 125 x 70mm 2 (•). 

number of branching events is strongly suppressed with 
decreasing d. This suppression, coupled with the disap- 
pearance of branch-lines with their arrival at the sample 
faces at Z = or d (see e.g. Fig. lc), causes the increas- 
ingly sporadic appearance of micro-branches as d — > AZ. 
Thus for d — > AZ (AZ ~ [im for PMMA and less for 
glass) micro-branching events are rare with crack dynam- 
ics becoming effectively two-dimensional. This suppres- 
sion of micro-branches explains how oscillations are so 
easily obtained for d ~ 0.15 — 0.2mm. 

Although rare, the micro-branching instability is ob- 
served in thin sheets both in straight and oscillatory 
propagation. When excited, micro-branches reduce the 
energy available to the main crack 2], and consequently 
reduce the main crack's velocity. If V falls below Vc, 
oscillations immediately cease, and the crack propa- 
gates along a straight trajectory (e.g. Fig. 5 in the 
2nd shaded region). When micro-branching is sporadic, 
micro-branches only momentarily reduce V to below Vc- 
This not only causes discontinuities in the oscillation 
phase but before and after a velocity drop we often ob- 
serve entirely different sinusoidal solutions (with signifi- 
cantly different wavelength and amplitude). When pro- 
longed micro-branching occurs, as characteristic of the 
thicker gels, V only momentarily surpasses Vc between 
micro-branching events. In this case, pronounced "stair- 
case" profiles are created by momentary diversions of the 
crack for less than a period followed by straight line tra- 
jectories. All of these effects are demonstrated in the 
single experiment, presented in Fig. 5, in which d is ta- 
pered as a function of X. The rarity of micro-branches in 
the initial stage (d < 0.35mm) enables rapid acceleration 
to V > Vc and the onset of oscillations. As the taper 
widens, micro-branching occurs; first sporadically (2nd 
shaded region) and later with prolonged branching (for 
d > 0.35mm), leading to stair-case profiles. 

In recent years, two oscillatory instabilities of cracks 
traveling in thin sheets have been reported. Both of these 
are fundamentally different from the oscillatory instabil- 
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theoretical description of the near vicinity of the crack 
tip. 
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FIG. 5: (color online) Tapering of d (0.19 < d < 0.43mm) as 
a function of X enables V > Vc for thicknesses (d > 0.35mm) 
normally dominated by micro-branches. (top) Schematic 
drawing of the XZ fracture surface indicating branch- line lo- 
cations and d(x). (center) Photograph of the crack profile 
along the XY plane with (bottom) the corresponding velocity 
profile. Shaded rectangles indicate regions of sporadic branch- 
lines (diamond shapes at top) and the resulting drops in V 
until their disappearance when the branch-line encounters the 
Z — or d planes. When V < Vc (Vc denoted by the dashed 
line) oscillations cease with no hysteresis and cracks continue 
along a straight-line trajectory (e.g. second shaded region). 
Prolonged branching (d > 0.35mm) yields "staircase" profiles. 



ity presented here. Deegan and coworkers have observed 
wavy cracks traveling at intersonic velocities in biaxially- 
stretched thin rubber sheets 14]. In our experiments, 
oscillations are observed at clearly subsonic velocities in 
pure uniaxial tension. Phase field models which do 
yield subsonic oscillations of crack trajectories, predict 
A to scale with the system size, also in sharp contrast 
with our results. Recent calculations ^(|, pairing LEFM 
with the Hodgdon and Sethna 3 path selection criterion, 
predict the onset of oscillatory instability at close to our 
measured values of Vc- In contrast to our measurements 
however, A is predicted to scale with the system size. 

Oscillatory cracks have also been observed when a rigid 
cutting tool is forced through a thin elastic sheet fl7j . 
These cracks travel at V « Cs with out-of-plane bend- 
ing playing an important role. In contrast, the instabil- 
ity presented here is in pure Mode I with no measurable 
out-of-plane deflection 18]. Our oscillations are triggered 
at Vc values that are near Cs, irrespective of the sam- 
ple's thickness. In addition, the crack dynamics in our 
experiments are consistent with the single-crack (Mode 
I) equation of motion predicted by LEFM until either 
Vc or the onset of micro-branching. We therefore con- 
clude that these observations indicate a previously unob- 
served instability in dynamic fracture whose character- 
istic time/length suggests that a new intrinsic scale is 
important to describe these dynamics. This time/length 
scale can not be explained in the framework of LEFM 
[l(| . underlining the necessity for a more fundamental 
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